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ABSTRACT 

While the satellites of the Milky Way (MW) have been shown to be largely con- 
sistent in terms of their mass contained within one half-light radius (M^aif) with 
a "universal" mass profile, a number of M31 satellites are found to be inconsistent 
with these relations, and seem kinematically colder in their central regions than their 
MW cousins. In this work, we present new kinematic and updated structural proper- 
ties for two M31 dSphs, And V and And VI using data from the Keck Low Resolu- 
tion Imaging Spectrograph (LRIS) and the DEep Imaging Multi-Object Spectrograph 
(DEIMOS) instruments and the Subaru Suprime-Cam imager. We measure systemic 
velocities of v r = —393.1 ± 4.2 km s -1 and —344.8 ± 2.5 kins -1 , and dispersions of 
a v = 11.5^4 4 kms -1 and a v — 9.4^2 4 kms -1 for And V and And VI respectively, 
meaning these two objects are consistent with the trends in <j v and rhaif set by their 
MW counterparts. We also investigate the nature of this scatter about the MW dSph 
mass profiles for the "Classical" (i.e. My < —8) MW and M31 dSphs. When com- 
paring both the "classical" MW and M3I dSphs to the best-fit mass profiles in the 
size-velocity dispersion plane, we find general scatter in both the positive (i.e. hotter) 
and negative (i.e. colder) directions from these profiles. However, barring one excep- 
tion (CVnl) only the M31 dSphs are found to scatter towards a colder regime, and, 
excepting the And I dSph, only MW objects scatter to hotter dispersions. The scatter 
for the combined population is greater than expected from measurement errors alone. 
We assess this divide in the context of the differing disc-to-halo mass (i.e. stars and 
baryons to total virial mass) ratios of the two hosts and argue that the underlying 
mass profiles for dSphs differ from galaxy to galaxy, and are modified by the baryonic 
component of the host. 



1 INTRODUCTION 


The past few years have been a revelation for the kine- 
matic pro perties of the dwarf s pheroidals (dSphs) of the Lo- 
cal Group. Istrigari et al.l (|2008l ) compiled kinematic data for 
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18 of the Milky Way dSph galaxies using a maximum like- 
lihood technique based on the Jeans equation to determine 
the masses for each of these systems within a 300 pc ra- 
dius (M3oo). This radius was chosen as the masses of these 
objects are best constrained within the region where there 
are tracers of the potential (i.e. stars) and for their sample, 
300 pc represented the average radius for this region. They 
determined that despite a luminous range of more than 4 or- 
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ders of magnitude, the objects were consistent with having a 
dynamical mass of 10 7 Mq within 300 pc of their centre, and 
declared this as a common mass scale for dSph galaxies. This 
characteristic mass scale h ad already b een observed for the 
brighter dSphs (My < -8. lMateolll99Sf ), but the consistency 
of the fainter objects was a surprise. Further, they showed 
that the dSphs were all consistent with having formed in 
stellar halos with total masses > 1O 9 M , which could im- 
plicate this as the cut-off mass for star formation within 
cold dark matter halos, or the minimum mass with which 
a dark m atter halo cou l d form . T his worked was th en ex- 
tended by I Walker et~ai1 (120091 1 and I Wolf et ail (|2010r ). both 
of whom show that the Milky Way (MW) dSph halos ex- 
hibit a correlation between global velocity dispersion and 
half-light radius, and as such, the central densities of their 
dark matter halos do not show signi ficant scatter o ver a 
large range of luminous scale-radii. I Walker et al.l (|2009l ) used 
this observation to postulate that the dSphs of the MW ap- 
pear to be consistent with having formed with a "Univer- 
sal" dar k matter halo mass profile. In addition, follow-up 
work by IWalker et ail (|2010h demonstrated that the halos 
of these dark matter dominated systems are also consis- 
te nt with the mean rotati on curve derived for spiral galaxies 
in lMcGaugh et al] {2007), indicating a constant dark mat- 
ter central surface density for galaxies ranging from MW 
dSphs to spirals. Recent results presented by both the Pan- 
Andromeda Archaeological Survey (PAndAS) and the Spec- 
troscopic and Photometric Landscape of Andromeda's Stel- 
lar Halo (SPLASH) groups have indicated significant dif- 
ferences between the kinematic temperatures (or velocity 
dispersions) of the dark matter halos of dSphs orbiting the 
MW vs. those orbiting Andromeda, with a number of the lat- 
ter being seemingly colder for a given h alf-light radius, and 
therefore less mas sive and less dense (|Collins et al.l l20ld ; 
iKalirai et aHl2010r i, signalling that this "Universal" mass 
profile may not extend to dSp hs found outside of the MW. 
They are also outliers to the iMcGaugh et al. (2007) rota- 
tion curve relations as discussed in lWalker et all (20101 This 
strongly suggests that while many dSphs appear similar in 
terms of their central densities, they are not all embedded 
within dark matter halos that follow a "Universal" profile. 

Interpreting the observed differences in the subhalos for 
these two populations is not trivial as it requires an under- 
standing of why some dSphs are born, or currently reside 
in, dark matter halos with differing physical and dynamical 
properties. One explanation could be that the underlying 
physical processes behind the formation of dSph halos dif- 
fers between the MW and M31; this clearly confronts our 
current cosmological paradigm for the growth of dark mat- 
ter halos. Another possibility is that the M31 subhalos could 
have formed later than the MW population when the Uni- 
verse was less dense, resulting in colder, less dense subhalos. 
The MW is thought to have had a more quiescent merger 
history than M31, accreting less baryons over cosmic time. 
This means that the MW could have 'formed' earlier (i.e. 
reached its half mass earlier) than M31. It is unclear how 
much later it would need to form to reproduce these re- 
sults so more modelling is require d to investigate this sce - 
nario fully. Finally, recent work by IPenarrubia et al.l (|2010T ) 
put forward a more physically motivated theory for this dif- 
ference. Using N-body simulations, they demonstrated that 
subhalos that evolve in host environments with disc-to-halo 



mass ratios of twice that measured in the MW naturally end 
up with less mass within their half-light radii, very similar to 
the observati ons of the M31 popula tion. This fits nicely with 
the results of iHammer et all (|2007T ). where they report that 
M31 has been more successful at accreting baryons than the 
MW, resulting in a stellar population that is ~ 2.5 times the 
mass of the MW stellar population. 

While there is a tendency for M31 dSphs to inhabit 
colder, lower mass halos, it is by no means observed through- 
out the population. Almost half of all spectroscopically sur- 
veyed M31 satellites (5 ou t of 12 - And I , And VII, And XIII 
And XV and And XVI iLetarte et ail 120091 ; ICollins et all 
2010; Kalirai et al.| [2010) are observed to have velocity dis- 
persions that are entirely consistent with their MW coun- 
terparts, suggesting an overall scatter within this popula- 
tion. Therefore, how significant is this subset of kinemati- 
cally colder M31 dSphs? Could this trend be the result of 
observational biases present in both systems? Due to our 
position in the MW, located within the disc at a distance 
of 785 kpc from M31, we are obviously not able to observe 
exactly the same subset of dSph galaxies in our own halo 
compared to that of M31. In the MW, we are hampered by 
obscuration from components of our own Galaxy that pre- 
vent an areal coverage of the system. However, for the region 
that has been covered in large all-sky surveys (such as SDSS) 
we are thought to be c omplete in our detections of bright, 
"clas sical" -type dSphs. l|Tollerud et al.l l2008; Ko posov et alj 
2008) In M31 the limiting factors in dwarf detection are 
attributed to the large distance between ourselves and our 
neighbour, which prevents us from observing objects with 
Mv > —6. We also struggle to identifying dwarf galaxies 
that lie close to the centre of M31 (at Rmzi < 40 kpc) as 
these objects are obscured (either directly or via projection 
effects) by the large disc component of the galaxy (traced 
out as far as 40 kpc kinematically and 25-30 kpc photometri- 
cally). We also observe more bright dSphs at large distances 
from the centre of M31 than we do in the MW, resulting in a 
mean distance of satellite from ho st of 184 kpc in M31 com- 
pared with 138 kpc in the MW (|McConnachie et all I20091 ; 
iRichardson et al.ll201lf ). 

The study of the intrinsic properties of dSphs has also 
been illuminated by the surprising behaviour that has been 
observed in some of the least luminous members of this pop- 
ulation, the "ultra-faint" dSphs. Since 2005, a large number 



of ve r y faint objects (—2 > My > —7, e.g. Willman et al.l 



20051; IBelokuroy et all 120061 . 120071 ; IZucker et all l2006bl lal: 



Belokurov et al 20081 ) have been discovered, and have re- 
vealed several unusual properties for the population. For 
example, they appear to be hugely dark matter dominated, 
resulting in a significant departure from the w ell-establishe d 
mass-luminosity relationship for such objects (lMateolll998l ). 
The exact morphological nature of these objects is also 
a subject of some controversy as they possess properties 
that are common to both dSphs and globular clusters and 
are therefore difficult to classify as eithe r, such as the 
unus u al objects, Willman 1 and Segue 1 (| Willman et all 
120051; IBelokuroy et all 120071; iNiederste-Ostholt et all \2WM , 
ISimon et al]|2010l ; IWillman et aill2010l '). They may also be 
far from dynamical equilibrium, and instead are undergo- 
ing extreme tidal disruption, inflating their observed veloc- 
ity dispersions a nd inferred mass-to-light ratios (e.g. Segue 
1 and Hercules, INiederste-Ostholt et all 120091 ; ISand et all 
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2009 and I Jin & Martini lioich . resulting in a misleading rep- 
resentation of the behaviour of these low-luminosity satel- 
lites. These faint objects are observed within a regime where 
inflation of the ve locity dispersion from binary star systems 
is no n-negligible |McConnachie fc Cote] |2010| ; iMinor et al.l 
l2010l ), meaning their dynamical dispersions may be lower 
than the values quoted in the literature, drastically altering 
the inferred M/L ratios. Therefore, if we restrict ourselves 
to comparing solely the dSphs whose kinematic properties 
are free from such uncertainties for the two host galaxies, 
does this variation persist? And what of its significance? 

To answer these questions, we have embarked upon 
a limited comparison of "classical" M31 and MW dSphs, 
where we define "classical" to be any object with Mv < 
— 7.9, with particular focus on two M31 dSphs that have re- 
cently been spectroscopically observed, And V and And VI, 
allowing us to present their kinematic properties for the first 
time. As And V and And VI represent the last two classical 
dSphs that belong to either the MW or M31 to have their ve- 
locity dispersions and masses measured, these observations 
make for a timely consideration of the whole ensemble of 
classical dSphs. 



2 OBSERVATIONS AND ANALYSIS 
TECHNIQUES 

2.1 And V 

2.1.1 Structural properties with Subaru Suprime-Cam 

Andromeda V is a bright M 31 satellite with Mv = —9.8 
|McConnachie fc Irwml2006ai ). an d was first discovered wit h 
the Second Palomar Sky Survey jArmandroff et all 1 19981 ). 
Its detailed structural properties we re assessed using Isaac 
Newton Telescope (INT) data by iMcConnachie fc Irwinl 
(2006a). This imaging was performed in Johnson V and 
Gunn i-band with the Wide Field Camera (WFC), and was 
deep enough to observe the top few magnitudes of the red gi- 
ant branch (RGB) for this object allowing the measurement 
of the surface brightness profile, intensity-weighted centre, 
position angle, ellipticity and scale-radii of the dSph using 
resolved counts of RGB stars. We summarize their results in 
Table[TJ During the nights of August 3-5 2005, deep Subaru 
Suprime-Cam imaging of a number of Andromeda dSphs, 
including And V were obtained (P.I. N. Arimoto) in the 
Cousins V— and I c — bands. Conditions were photometric 
throughout with typical seeing of 0.5". Full details regard- 
ing the observing strat egy and reduction technique s for this 
survey are outlined in IMcConnachie et al.l (|2007l ). but we 
briefly summarize here. Objects were typically observed in 
5 x 440 and 20 x 240 second exposures in V and I c bands re- 
spectively, allowing equivalent depths to be reached in both 
bands. Data were processed using the CASU photometric 
pipeline llrwin fc Lewis! l|200lf) . which de-biased, flat-fielded, 
trimmed and gain corrected the images. Then, a catalogue 
was generated for each image frame, morphologically clas- 
sifying each object as stellar, non-stellar and noise-like. We 
show the colour magnitude diagram (CMD) generated from 
the stellar-classified objects from these observations, where 
we have corrected for e xtinc tion and reddening using the 
maps of ISchlegel et~ai1 (|l998T ). in Fig. □ We now use this 
dataset to revise the structural properties of And V. To 



do this, we follow the same s tar count method detailed in 
IMcConnachie fc Irwinl (|2006bh so that we can make a like- 
for-like comparison of the Subaru and INT data. Briefly, 
we construct an isopleth map of a 60' x 60' field from the 
Subaru data centred on And V, then determine the cen- 
tre of gravity, ao,So, position angle 6 = 32 it 2° (measured 
from East to North), and ellipticity, e = 0.17 ± 0.02, for 
each isophote by using the int ensity-weighted moments and 
IMcConnachie fc Irw in (2006b) equations 2 and 3. We then 
construct a background corrected radial profile for And V, 
where a background level is estimated by measuring the 
average number of stars per arcmin within circular annuli 
located beyond the tidal rad ius of And V (as derived in 
IMcConnachie fc Irwinl l2006al ). We then subtract this aver- 
age background from star counts performed in elliptical an- 
nuli based on the centre of the dSph, and fit the result- 
ing profile with exponential, Plummer and King profiles. 
Both the radial profile and best-fit models are displayed 
in Fig. [2] From this, we estimate the half-light radius (us- 
ing our result for the exponential scale-radius, r e = 0.88') 
to be rhaif = 1-3' and (for the King model) tidal radius 
of r t = 5.4' Using the dis tance modulus of 24.44 (774 
kpc lMcConnachie et akll2005l ) this gives r ha i f = 292 ±22 pc 
and r t = 1.2 ± 0.1 kpc. Our results are also summarised in 
Table[TJ 

We inspected the photometric metallicities of the And 
V stars using the Dartmouth isochrone models (Dotter et 
al. 2008). We select an age of 10 Gyrs and [a/Fe]=+0.2 as 
numerous studies of dwarf spheroidals have shown them to 
be composed of old stellar populations that are enhanced in 
a-elements. We interpolate between these isochrones on a 
fine grid, and present the resulting metallicity distribution 
function (MDF) in Fig. [3] We find a mean metallicity of 
[Fe/H]=-1.6 for the And V members and a dispersion of 0.3 
dex. However, as our imaging is not deep enough to realise 
the main sequence turn-off (MSTO) of And V, we cannot 
reliably ascertain the precise ages and a-abundances of the 
stellar populations. The effect of increasing (decreasing) the 
age used for our isochrones by ~ 2Gyr results in a shift of 
—0.1 (+0.1) dex to our mean metallicities. Similarly, and 
increase (decrease) in our assumed value of [a/Fe] by 0.2 
dex shifts our mean metallicity by —0.1 (+0.1) dex. This 
results in a combined uncertainty of ±0.2 dex. 

2.1.2 Kinematic properties with Keck I LRIS 

And V was observed using the Low Resolution Imaging 
Spectrograph (LRIS), situated on the Cassegrain focus of 
the Keck I telescope on Mauna Kea, on 16th August 2009. 
The 831/8200 grating was employed with the /-band fil- 
ter on the LRIS Red II detector system, giving a resolution 

o 

of ~ 3.0 A (R~ 2800). Our observations were taken in an 
average of 0.7" seeing, covering a wavelength range of 6900- 

o 

9000 A, in the region of the Calcium Triplet (Ca II). We took 
4x15 minute exposure, which resulted in typical signal-to- 

o 

noise ratios of S:N=3-10 A" 1 . Targets were selected using 
the Subaru CMD shown in Fig. [TJ prioritising stars with 
20.3 < i < 23.0 from within the locus of the clearly defined 
the RGB. In total, we observed 45 science targets within the 
And V mask, 38 of which reduce successfully. 

To reduce the LRIS data, we used the IRAF 
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Figure 1. Extinction corrected CMDs of And V and And VI (Subaru, V and I c ) for all stellar objects within 1 X r^if of the dSph 
centres (small black points). In both cases, the RGB is clearly seen. Red circles indicate confirmed spectroscopic members for each dSph 
with S:N> 3, blue circles indicate confirmed spectroscopic members for each dSph with S:N< 3 and cyan circles indicate stars with 
velocities that are consistent with the dSph systemic velocity, but are potential M31 halo contaminants. Larger black points indicate 
stars that were observed spectroscopically, but whose velocities are inconsistent with the systemic velocities for And V and And VI. 
Error bars show the average lc uncertainties in the photometry at each magnitude level 



noao.twodspec and noao.onedspec packages. The spectra 
were flat-fielded, de-biased and wavelength calibrated, as 
well as cleaned of skylines and cosmic rays. We then per- 
formed heliocentric velocity corrections to each of our ob- 
served stars. Our velocities and errors were derived from 

o 

the Ca II triplet lines locat ed at ~ 850 A, using the 
same technique described in ICollins et ail (|201Cj ). Briefly, 
we used an error-weighted cross-correlation technique with 
a model template of the Ca II feature at the rest-frame 
wavelength positions of the triplet lines. We repeated the 
cross-correlation 1000 times, adding random Poisson noise 
to our spectra on each occasion, and took the average and 
standard deviations of these as our velocities and associated 
errors. We then combine these errors with ones derived by 
performing separate cross-correlations to each of the Ca II 
lines in turn, using the dispersion of the resulting velocities. 
This results in typical errors of 5 — 15kms _1 , with a mean 
for our sample of 6.8 km s -1 . 

2.2 And VI 

2.2.1 Structural properties with Subaru Suprime-Cam 

And VI is one of the more lumin ous of the M31 clas- 
sical dSphs, with My — —11.5 (|McConnachie fc Irwin! 
l2006al ), and it too was discovered using t he Second Palo- 
mar Sky Survey l|Armandroff et ail [l999) . Its structural 
properties were assessed from INT WFC photometry by 
iMcConnachie fc Irwinl (|2006bh in the same manner as And 
V. These results are summarised in Tabled And VI was also 
observed in our Subaru Suprime-Cam survey, and the data 
were taken and reduced in the same way as described above. 
We display the CMD of stellar objects from this data within 
1 x rhai f of And VI in Fig. [1] The data start to become in- 
complete at a S:N of 10, which corresponds to V = / ~ 25.5. 



In this section, we re-derive the structural parameters 
for And VI, using this Subaru data and the same approach as 
detailed for And V. Again, from the star counts and isopleths 
maps we find a good agreement with those of the INT study. 
Our results are summarised in Table(2] We calculate ruaij = 
1.9' (from our exponent ial scale radius), rt = 7.0' , e = 0.39 
and 6 = 164° . Using the lMcConnachie et all (|2005l ) distance 
for And VI of 783 kpc (distance modulus of 24.47) we obtain 
physical scale radii of rhaif = 440 ± 16 pc and rt = 1.6 ± 
0.2 kpc. 

We perform the same isochrone analysis for And VI as 
that described in §2.1.1, and the resulting MDF is displayed 
in Fig. [3] The same difficulties in ascertaining robust photo- 
metric metallicities for And V are experienced here, as the 
Subaru data do not reach the MSTO. We find a median 
metallicity of [Fe/H]= —1.1 ± 0.3, where the error repre- 
sents the dispersion of the population. This dispersion is 
the same as we measure for And V, despite And VI having 
a much broader RGB. And V is more metal poor than And 
VI with an average metallicity of [Fe/H] = —1.6. As you 
decrease in metallicity, isochrones for a given age and en- 
abundance begin to bunch together in colour space (V — I), 
meaning smaller differences in colour translate to larger dif- 
ferential metallicities, driving a larger metallicity dispersion 
for a smaller colour range. 

2.2.2 Kinematic properties with Keck II DEIMOS 

The DEep-Imaging Multi-Object Spectrograph (DEIMOS), 
situated on the Nasmyth focus of the Keck II telescope is an 
ideal instrument for obtaining medium resolution (R~ 6000) 
spectra of multiple, faint stellar targets in the M31 dSphs. 
The data for And VI were taken on the nights of 17 -19 
September 2009 in photometric conditions and with < 1" 
seeing. Our chosen instrumental setting covered a wave- 
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Figure 2. Radial profile for And V (left) and And VI (right) constructed from aperture star counts in the region of both dSphs. In each 
case, an average background number density was estimated from an annulus located between 8' — 12' from the centre of each object, 
and the results were subtracted from our star counts. We overlay the best fit Exponential (blue short dashed line), Plummer (magenta 
long-dash line) and King (green solid line) profiles in each case also. 
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Figure 3. The MDF for And V and And VI, constructed from metallicities f or stars within 2rh a i f of the centre of each dSph. These 
were calculated by interpolating between a fine grid of Dartmouth isochrones llDotter et al.|[2~00Sl) with [ct/Fe]=+0.2 and an age of 10 
Gyr. The mean metallicity is found to be [Fe/H] = —1.6 ± 0.3 and [Fe/H] = —1.1 ± 0.3 dex for And V and And VI respectively where 
the error represents the dispersion of the population. 



length range of 5600-9800 A and for our exposures we 
implemented 3x30 minute integrations, and employed the 
1200 line/mm grating, giving us a spectral resolution of 

o 

~ 1.37 A. The spectra from this setup typically possess S:N 

o 

of > 3 A - . Due to the position of And VI in the sky, this 
object was observed at the beginning of each night meaning 
that it suffered from a high airmass of 2.2, so the S:N for 
these stars are lower than expected for comparable length 
exposures taken at a lower airmass. The typical errors within 



the mask range from 5 - 
pie of 6.2 kms~ . 



12 km s , with a mean for our sam- 



2.3 Determining membership 

Before we can analyse the kinematic properties of these ob- 
jects, it is important to correctly determine the bona-fide 
members of the dSphs, and to eliminate stars which belong 
to the MW foreground (mostly v r > — 160 km s -1 ) or the 
M31 halo (v r ~ -300 kms" 1 ). As both And V and And VI 
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Table 1. Structural properties of And V derived in McConnachie 
et al. 2006a and this work. 



Table 2. Structural properties of And VI derived in McConnachie 
et al. 2006a and this work. 



Property 


M06a 


This work 


Property 


M06a 


This work 


ctO (h:m:s) 


01:10:17.0 


01:10:17.9 


ao (h:m:s) 


23:51:46.9 


23:51:47.3 


So (° :':") 


+47:37:46 


+47:37:38.0 


So (° :':") 


+24:34:57 


+24:34:52.5 


e (°) 


32±3 


32+2 


(°) 


163+3 


164+2 


e 


0.18±0.03 


0.17+0.02 


V 


0.41+0.03 


0.39+0.02 


r e (arcmin) 


0.86±0.05 


0.88+0.09 


r e (arcmin) 


1.2+0.04 


1.1+0.04 


r p (arcmin) 


1.56±0.08 


1.59+0.1 


r p (arcmin) 


2.15+0.08 


2.3+0.08 


r c (arcmin) 


1.2+0.2 


1.2+0.1 


r c (arcmin) 


2.1+0.2 


1.5+0.1 


rt (arcmin) 


5.3+0.4 


5.4+0.8 


rt (arcmin) 


6.2+0.4 


7.0+0.8 


rt (kpc) 


1.2+0.2 


1.2+0.2 


rt (kpc) 


1.4+0.1 


1.6+0.2 


rhalS (arcmin) 


1.3 


1.3+0.1 


rhais (arcmin) 


1.8 


1.9+0.07 


Thalf (P c ) 


300 


292+22 


ThalS (P c ) 


420 


440+16 




-1.6 


-1.6+0.3 




-1.3 


-1.1+0.3 



(o) E rror represents dispersion. 

lie at a large projected ra dius from the centre of M31 (~ 110 
and 270 kpc respectively, lArmandroff et al.| [l998. 1999), the 
density of stars, and thus the number of contaminants, be- 
longing to the M31 halo should be very low for these ob- 
jects, but should still be considered. As And V is located at 
a higher latitude than And VI, this dwarf is more likely to 
experience contamination from the MW, and so this must 
be treated carefully.In order to minimise the contamination 
to our sample, we require the following conditions to be met 
by bona-fide members: 

• Stars must fall on the RGB of dSph, as defined by colour 
cuts in V and /- bands for And V and And VI. 

• Stars must sit within 2rhais of the dSph centre. Stars at 
larger radii that satisfy all other criteria barring this are 
classified as 'tentative members', and their properties and 
probabilities of being members will be discussed. 

• Stars must show low Na I doublet absorptions (i.e. 
EWjva/ < 1-4). Strong Na I absorption may indicate that 
the observed star is a foreground dwarf star, not an M31 
RGB star. 

• After deriving a systemic velocity, any star within the 
sample with a velocity that is > 3<r from the systemic will 
not be considered a member. The kinematic properties of the 
dSph will then be re-derived iteratively in this manner until 
only stars with velocities within 3er of the derived systemic 
velocity remain. 



3 THE KINEMATICS OF And V AND VI 

3.1 Systemic velocities and velocity dispersions 

Using the criterion laid out in § 2.3 we identify the most- 
probable member stars for each of the dSphs observed with 
LRIS and DEIMOS. In Figure Q] we display the kinematic 
properties of all the stars observed in each slit mask. The 
top panel in each case shows a histogram of the velocities of 
observed stars. The heavy red histogram highlights the kine- 
matic location of each dSph, which are immediately obvious 
as cold over-densities of stars in velocity space. The cen- 
tral panel shows velocity vs. distance from the dwarf centre, 
where the 1, 2, 3 and 4x Thais are marked as dashed lines. 
The lower panel shows the velocity as a function of photo- 
metric metallicity. Not all stars within the mask have 



(a) j5 rror represents dispersion. 

photometric metallicities as their colours place them 
outside the parameter space covered by our chosen 
isochrones. Such drop outs are likely Galactic con- 
taminants and not true M31 RGB stars. Cutting on 
velocity, distance, Nal absorption and metallicity (proxy for 
position on RGB), we identify 14 and 26 secure members 
per dSph for And V and And VI respectively, as well as 
2 and 7 (4 of which have S:N> 3) tentative members be- 
yond 2 x Thais that fall within the cold peak, although we 
treat these cautiously as at the systemic velocities for these 
objects (~ — 400 kms -1 and ~ — 350 kms -1 ) fall within the 
regime of the M31 halo. In the case of And V, two stars found 
within the velocity peak do not fall on the RGB of the dwarf 
(shown as cyan points in Fig.[T]), seemingly more metal poor 
than the And V population, making them likely halo inter- 
lopers. For And VI we also identify 2 outliers from the CMD 
which again are likely halo contaminants. In both cases, the 
outliers show no significant Na I absorption, strengthening 
our assumption that these are M31 halo RGB stars, rather 
than Galactic halo dwarfs. We show the kinematic proper- 
ties of all confirmed and tentative members, as well as halo 
contaminants for And V and And VI in Tables [3] and [4] re- 
spectively. 

With secure candidates in hand, we now determine 
the systemic velocities (v r ) and velocity dispersions (a v ) 
for each dSp h using the maximum likelihood approach of 
iMartin et al.l (2007). This is the s ame technique as we 
have used in our previous work (e.g. Chapman et al.l 120061 : 



iLetarte et al.|[2003 : ICollins et al.ll201C ). and we briefly sum- 
marise here. We calculate v r and a by sampling a coarse 
grid in (v r ,a) space and determining the parameter values 
that maximise the likelihood function (ML), denned as: 



ML(v r 



exp 



O"tot 



(1) 



with N the number of stars in the sample, « r ,i the radial 
velocity measured for the i th secure member star, v err ^ 



the corresponding uncertainty and a to t = v/ 'a 2 + v 2 r r,i- ^ n 
this way, we are able to separate the intrinsic dispersion of 
the dSph from the dispersion introduced by our measure- 
ment uncertainties. We display the one dimensional likeli- 
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Figure 4. The velocities, distances from centre and photometric metallicities of all stars within the LRIS and DEIMOS fields. Some stars 
do not have photometric metallicities as their colours place them outside the range of isochrones used. Red points indicate likely dSph 
members and the red dashed lines represent 1, 2, 3 and 4 times the r/j of each dSph. And V and And VI appear as cold over-densities 
of stars at ~ — 400 kms -1 and ~ — 350 kms - 1 respectively. 



Table 3. Kinematic properties of observed members of And V 

Confirmed members 



«20oo(hh:nrm:ss) 


feooo 


V 


1 


v r ( km s 1 ) 


S:N ( A -1 ) 




[Fe/H] spec 


1:10:18.62 


47°38'9.4" 


22.544 


21.405 


-365.6±12.6 


6.5 


-1.65 


-2.1±0.3 


1:10:12.19 


47°37'30.6" 


21.95 


20.581 


-377.3±21 


5.1 


-1.55 


-1.3±0.3 


1:10:20.65 


47°38'15.7" 


21.979 


20.735 


-408. 2±4 


4.6 


-1.85 


-2.2±0.4 


1:10:20.03 


47°37'10.1" 


22.349 


21.101 


-400.4±12.8 


4.4 


-1.48 


-1.2±0.3 


1:10:16.32 


47°37'37.4" 


22.803 


21.815 


-396.2±2.21 


4.3 


-2.15 


-1.6±0.7 


1:10:22.22 


47°38'52.1" 


22.685 


21.571 


-405.7±11.5 


4.2 


-1.55 


-2.1±0.3 


1:10:25.14 


47°38'9.8" 


22.875 


21.873 


-401.6±1.94 


3.9 


-2.05 


-2.0±0.3 


1:10:17.94 


47°37'16.2" 


22.05 


20.804 


-366.4±11.7 


3.5 


-1.78 


-1.0±0.4 


1:10:06.71 


47°35'26.9" 


22.218 


20.986 


-423.8±1.28 


3.4 


-1.65 


-1.8±0.4 


1:10:09.04 


47°37'20.8" 


22.294 


21.125 


-388.1±4.1 


3.4 


-1.78 


-1.6±0.5 


1:10:24.17 


47°37'47.1" 


21.323 


19.909 


-408.3±9.04 


2.8 


-2.12 


-1.7±0.7 


1:10:21.56 


47°37'25.8" 


22.152 


20.86 


-382.3±6.76 


2.2 


-1.55 


-1.8±0.5 


1:10:05.54 


47°36'41.6" 


22.178 


20.871 


-406.5±3.6 


2.0 


-1.48 


-1.4±0.5 


1:10:02.38 


47°37'48.5" 


22.84 


21.808 


-371.3±5.08 


1.7 


-1.85 


-2.2±0.5 




Tentative members ( 


.e. stars at > 2 X Thalf) 






O200o(hh:mni:ss) 


feooo 


V 


I 


v r ( km s 1 ) 


S:N ( A -1 ) 


Fe/H]$ ot 


[Fe/H] spec 


1:10:07.45 


47°35'48.0" 


21.845 


20.526 


-371.2±2.34 


3.9 


-1.9 


-1.9±0.4 


1:10:38.56 


47°38'22.1" 


22.51 


21.346 


-409.0±4.38 


3.0 


-1.55 


-2.0±0.4 








Outliers 








O200o(hh:rrim:ss) 


feooo 


V 


1 


v r ( km s 1 ) 


S:N ( A" 1 ) 




[F^/H] spec 


1:10:19.23 


47°39'4.1" 


21.748 


20.563 


-390.2±1.32 


3.8 


-2.42 


-1.3±0.5 


1:10:13.94 


47°37'9.0 " 


22.067 


20.989 


-381.3±2.03 


3.7 


-2.48 


-1.6±0.7 



(a) j) er i vec i from lDotter et al] j2008! ) isochrones with [a/Fe]=+0.2, age=10 Gyrs. Typical errors of ±0.2 dex. 



hood distributions for v r and a v in Fig. [5] where the dashed 
lines represent the conventional 1, 2 and 3cr (68%, 95% and 
99.7%) uncertainties on the values. In the case of the And 
VI data, we have only included the 12 member stars with 

o 

S:N> 3 A in this analysis to avoid introducing large uncer- 



tainties into their calculated properties. We determine sys- 
temic velocities of v r = -393.1±4.2 and -344.8±3.4kms _1 
and velocity dispersions of a v — 11. 5l|' 4 and 9. 41^4 kms -1 
for And V and VI respectively. The velocities of these two 
galaxies (but not dispersions) were previously reported in 



8 M. L. M. Collins et al. 

Table 4. Kinematic properties of observed members of And VI 

Confirmed members 



Q 2000 (nn:mm:ss) 


O2000 


i / 

V 


I 


/ l —W 

v r ( km s ) 


o 

Q.TVT / A — 1 \ 

SUN 1 A ) 


rt.vwtn( a ) 


fey, /XT] 


23:51:44.45 


24°37'48.80" 


22.12 


20.96 


-347.25±1.73 


5.02 


-1.22 


-2.04 ±0.8 


23:51:50.74 


24°33 '57.70" 


23.76 


22.65 


-340.84±2.36 


4.57 


-0.70 


-0.96±0.5 


23:51:45.55 


24° 36' 3.80" 


22.15 


21.09 


-324.21±2.57 


4.38 


-1.35 


-0.71±0.2 


23:51:52.41 


24°33'34.70" 


22.83 


21.63 


-340.42±2.14 


4.04 


-0.80 


-1.41±0.5 


23:51:51.52 


24°36'23.50" 


22.46 


21.29 


-345.70±3.48 


3.84 


-1.04 


-0.98±0.2 


23:51:46.94 


24°35'40.70" 


22.33 


21.19 


-326.74±3.01 


3.81 


-1.12 


-1.71±0.4 


23:51:48.19 


24°34'27.10" 


24.98 


24.09 


-337.44±2.64 


3.57 


-0.92 


-1.15±0.4 


23:51:41.15 


24°36' 33.00" 


22.38 


21.29 


-341.16±3.32 


3.57 


-1.15 


-0.78±0.3 


23:51:46.60 


24°37'11.90" 


22.35 


21.27 


-323.01±2.88 


3.56 


-1.18 


-2.45±2.2 


23:51:50.85 


24°36'18.50" 


22.46 


21.18 


-336.70±2.58 


3.55 


-0.99 


-1.81 ±0.6 


23:51:44.38 


24° 34' 16. 10" 


24.06 


23.28 


-350.49±3.40 


3.20 


-0.73 


-0.78 ±0.4 


23:51:46.21 


24°37'37.10" 


23.61 


22.91 


-347.09±5.88 


3.10 


-1.18 


-1.37±0.3 


23:51:46.19 


24°35'54.60" 


25.70 


24.90 


-332.98±6.07 


2.08 


-1.15 


-1.09±0.4 


23:51:41.12 


24°37' 43.50" 


23.00 


22.09 


-308.26±8.37 


2.00 


-1.05 


-1.59 ±3.2 


23-51 -54 77 


24°28'53.50" 


23.19 


22.20 


-339 81 +6 66 


1.90 


-0.80 


-1.87±0.9 


23:51:46.91 


24°33'53.00" 


23.69 


22.70 


-314.51±6.52 


1.89 


-0.95 


-0.95±0.7 


23:51:54.79 


24°35'50.30" 


22.88 


21.96 


-328.10±7.15 


1.86 


-1.10 


-1.41 ±1.1 


23:51:44.42 


24°37'26.30" 


23.38 


22.46 


-348.25±6.95 


1.46 


-0.83 


1.57± 11.1 


23:51:54.23 


24°36'42.80" 


22.80 


21.94 


-351.77±12.16 


1.43 


-1.25 


-1.7±0.9 


23:51:44.13 


24°35'45.10" 


23.37 


22.44 


-332.82±8.19 


1.39 


-0.80 


-1.14±0.6 


23:51:49.63 


24°32'42.20" 


22.98 


22.12 


-354.98±8.92 


1.35 


-1.15 


-2.16±11.4 


23:51:44.66 


24°36'37.90" 


24.83 


23.91 


-340.10±35.83 


1.21 


-1.15 


-0.54±0.5 


23:51:46.70 


24°35'16.50" 


24.09 


23.49 


-354.44±11.69 


1.33 


-1.25 


-1.88±0.7 


23:51:45.61 


24°36'47.10" 


22.92 


21.96 


-362.91±8.38 


1.11 


-1.02 


-1.16±0.9 


23:51:44.69 


24°34' 42.30" 


25.44 


24.59 


-384.50±9.38 


0.89 


-0.90 


-0.87±0.9 


23:51:41.55 


24°36'28.10" 


23.27 


22.44 


-320.62±14.57 


0.87 


-1.05 


-1.03±0.9 




Tentative members ( 


i.e. stars at > 2 


X r ha lf) 






"200o(hh:mm:ss) 


<52000 


V 


1 


» r (km s — 1 ) 


O 

S:N ( A" 1 ) 




Fe/H] spec 


23:51:47.00 


24°32'56.00" 


23.56 


22.62 


-363.18±1.85 


4.71 


-0.68 


-1.20±0.3 


23:51:45.80 


24°32'31.20" 


22.76 


21.71 


-326.89±5.04 


3.29 


-1 


-1.41± 0.4 


23:51:54.12 


24°31'36.80" 


23.30 


22.34 


-348.75±5.65 


3.75 


-0.78 


-1.36±0.5 


23:51:49.59 


24°31'27.20" 


23.38 


22.23 


-345.96±8.11 


3.51 


-0.57 


-1.10±0.7 


23:51:49.30 


24°32'26.80" 


23.16 


22.11 


-343.06±5.88 


1.46 


-0.74 


-2.06±0.7 


23:51:45.30 


24°30'18.50" 


23.34 


22.24 


-351.28± 8.70 


1.15 


-0.60 


-1.34±0.8 


23:51:47.81 


24°31'18.10" 


23.29 


22.29 


-320.05±11.83 


1.10 


-0.73 


-1.97±0.5 








Outliers 








O200o(hh:mrn:ss) 


#2000 


V 


1 


v r ( km s 1 ) 


S:N ( A" 1 ) 




[Fe/H] spec 


23:51:44.87 


24°37'31.90" 


21.70 


20.62 


-339.94±1.55 


6.39 


NA 


-1.65±0.3 


23:51:44.69 


24°34'42.30" 


24.19 


23.06 


-384.50±9.38 


0.89 


-0.20 


-1.20±1.3 



(a) Derived from lDotter et all i2008l ) isochrones with [a/Fe]=±0.2, age=10 Gyrs. Typical errors of ±0.2 dex. 



iGuhathakurta et al.1 (|200Cj ) to be v r = -387.0 ± 4.0 km s" 1 
and v r = -340.7±2.9kms _1 . which a eree with our findings 
within the associated la errors. We summarise our results 
in Tabled 



3.1.1 Masses and dark matter content 

As dSph galaxies are dispersion supported objects, we can 
use their internal velocity dispersions to measure a mass 
for the system, and infer how dark matter dominated they 
are. There are severa l methods in the literature fo r this 
(e.g. Illlingworthl Il976r . iRichstone fc Tremainel Il98rj ). how- 
ever these methods make the assumption that mass follows 
light, something we know to be incorrect from measurements 



Table 5. Kinematic properties of And V and And VI. 



Property 



And V 



And VI 



Vr ( km s 1 ) 


-393.1±4.2 


-344.8±3.4 


(j„(kms _1 ) 


Li-5±J;1 


g^l^kms- 1 


[Fe/H] spec 


-1.6±0.3 


-1.3 ±0.14 


M ha lf (M ) 


2.3+J 3 x 107 


2-i±S:g x io 7 


[M/L] kalf 


78ltJ 


in o+5.9 
iZ -' r> -4.7 



of many dSph systems (e.g. IWalker et al ] |2007l . l2009l . whose 
high masses cannot b e explained by the luminous matter 
only. Recent work by IWalker et al.l (|2009h has shown that 
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Figure 5. ML distributions for the systemic velocities and velocity dispersions of And V (left) and And VI (right). Dashed lines represent 
the 1, 2 and 3a uncertainties on the peak values. We measure v r = —393.1 ±4. 2 km s — 1 and — 344.8±2.5 km s — 1 , ando\, = 11. 5jV 4 km s~ 1 
and a v = 9.4^ 4 km s _1 for And V and And VI respectively. 



the mass contained within the half-light radius (Mhaif) of 
these objects can be reliably estimated using the following 
formula: 



Mhalf = fir half&v, half 

where fj, = 580 M kpc _1 krn _1 : 



(2) 



~s~, Thai; is the half-light ra- 
dius in kpc and cr Vt h a lf is the velocity dispersion within the 
half-light radius. This simple estimator assumes that the 
stellar component is distributed as a Plummer sphere with 
an isotropic velocity distribution and a velocity dispersion 
remains constant throughout the system. If we apply this es- 
timator to And V and And VI we obtain the following values 
for M half : 2.3tH x 1q7 M o for And V and 2.1+J-jj x 10 7 M 
for And VI, also summarized in Table [5] From this, it is 
trivial to estimate the mass-to-light ratios for these objects. 
The V— band luminosities for And V and And VI are Ly — 
5.75 x 10 6 L and L v = 3.40 x 1O 7 L , and we obtain the 
following values: \M/L] h aif =7 8± 5 A and [M / 'L] half =12. 3t 4 ;? 
(IMcConnachie fc Irwin 2006a). This demonstrates that each 



of our objects is likely to be dark matter dominated. 



3.1.2 Metallicities 

In § 2, we used Subaru photometry to measure the photo- 
metric m etallicities of the R GB stars in And V and And 
VI. Using iDotter et alJ i|2008r i isochrones with [a/Fe]=±0.2 
and an age of 10 Gyrs, we deduced average metallicities of 
[Fe/H]= -1.8 and [Fe/H]= -1.1 for And V and And VI 
respectively. While these values give us a good sense of the 
relative metallicities of these objects, our photometric data 
do not go deep enough to realise the MSTO of these dSphs, 
preventing us from deducing accurate ages for the stellar 
populations within. We are also unable to reliably discern 
their a-abundances, and these factors leave us exposed to 
the age-[a/Fe]-metallicity degeneracy, which attributes an 
error on our estimates of ±0.2 dex. To gain further insight 
into the average metallicities of these objects, we can turn 
to the spectra of the members of each dSph. 

As discussed in § 2, for both And V and And VI, 
our observational setup was such that we observed our 



stars in the wavelength regime of the Ca II triplet. The 
equivalent widths of this strong abs orption feature have 
been shown by numerou s authors (e.g. iBattaglia et al.ll2008l : 
IStarkenburg et all [2OI0I ) to be a very good proxy for the 
mea surement of [Fe/H], dow n to values as low as [Fe/H] = 
-4 jStarkenburg et all 120101 ) . Thus, we are able to use this 
feature to measure the metallicities for the individual stars 
of And V and And VI as well as for a composite spectrum 
for each, where we perform an error weighted co-addition 
of the spectra of the confirmed members of each satellite. 
To calculate metallicities from the equivalent widths of the 
3 lines (Ca IIg4gg,Ca Il8542,Ca IIg662) we use the following 
relation: 

[Fe/H] = -2.66 + QA2[EW - 0.64{V HB - Vrgb)} (3) 

where EW = 0.5 Ca II 8 498+Ca IIs542+0.6 Ca II8662, Vhb 
is the F-band magnitude o f the HB at the distance of M31 
(25.17. lHarbeck et al.ll2005l ) and Vrgb is the l/-band mag- 
nitude of the star (or weighted average V-band magnitude 
of the sample for the composite) . 

In Fig.[6]we show a number of spectra for our And V and 
And VI members stars spanning a range of luminosities and 
S:N, along with their calculated spectroscopic metallicities. 
When calculating errors we not only include the errors on the 
fit itself, but systematic errors in measuring the continuum 
level and the relationship between the equivalent widths and 
the measurement of [Fe/H]. As can be seen, for the star- by- 
star metallicities our [Fe/H] measurements carry significant 
errors (of order 0.3-0.5 dex or greater), as precision measure- 
ments of equivalent widths and the continuum of each star 
are difficult to accurately assess in this low S:N regime. This 
inaccuracy can be seen in Tables l3l and l4l where the 
spectroscopic metallicities are often very different 
from the photometrically derived metallicities which 
typically have errors of only ±0.2 dex. From measuring 
the metallicities of the individual members, we find a range 
of metallicities for And V from the spectra of [Fe/H]= —2.2 
- [Fe/H]= -1.2, with a mean of [Fe/H]= -1.6. For And 
VI we find a range of [Fe/H]= -0.9 - [Fe/H]= -1.8 with a 
mean of [Fe/H]= —1.4. With the co-added spectra, in Fig. [7] 
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And V 



And VI 




408.2 km/s [Fe/H] = -1.8±0.3 

.J I I I I I I I I I I I I 




if =-098.2 km/si [Fe/H] = -1.5±0.4 

i I i i , J_J i i , , I , i i i 

6500 8550 8600 6650 

Wavelength (A) 




[Fe/^] = -1.3j:0f3 
8550 8600 8650 

Wavelength (A) 



Figure 6. 5 spectra from individual And V (left) and And VI (right), shifted to the rest-frame and normalised in flux. These span 
a range of S:N to demonstrate the quality of our datasets. We indicate the positions of the Ca II features with red dashed line. The 
velocities and spectroscopic metallicities of each star are shown also. 



0.0 I L_ 

8450 



And V 

[Fe/H] = -1.6 ± 0.28 




8500 8550 8600 

Wavelength (Angstroms) 




E 

n 0.6 



0.4 

0.2 



0.0 I L_ 

8450 



8500 8550 8600 

Wavelength (Angstroms) 



Figure 7. Co-added spectra for And V (left) and And VI (right), constructed from member stars of each with S:N> 3. We use the Ca 
II triplet lines to estimate the average metallicity for each dSph and find [Fe/H] = —1.6 ± 0.3 and [Fe/H] = —1.3 ± 0.14 for And V and 
And VI respectively. The larger errors in our And V estimates are due to the lower resolution of our LRIS data and the fact that we 
exclude the 3rd line of the triplet from our calculation of [Fe/H] as it has been artificially broadened by skylines. 



we derive metallicities of [Fe/H] = —1.6 ± 0.3 for And V 
and [Fe/H] = -1.3 ± 0.15 for And VI. The larger errors on 
the [Fe/H] for And V are the result of a number of factors. 
First, we have a lower number of member stars for this dSph, 
meaning the noise in the resulting spectrum is greater. Sec- 
ond, the spectroscopic resolution of our LRIS setup is lower 
than that of our DEIMOS setup, which can cause broaden- 
ing effects in the Ca II lines. Finally, as can be seen in both 
the individual and combined spectra for this object, the 3rd 
line of the Ca II triplet has been broadened significantly by 
OH absorption lines in this regime, and so we have excluded 
this line from our metallicity calculations, using only the 
first and second lines. This is done with the knowledge that 



the ratio of equivalent widths of the Ca II lines when they 
are n ot saturated should be 0.4:1:0.75 IjStarkenburg et all 
2010), so we can alter the ratios of the Ca II8498 lines used 
in eqn. 3 from 0.5 to 1.5. We also see broadening in a num- 
ber of the And VI spectra (as can be seen in Fig [5}, but it 
is less prevalent throughout the sample. In the case of both 
And V and And VI, these spectroscopic values agree very 
well with those we derived from the photometry. They also 
conform to the observed trend of decreasing metallicity with 
decreasing lu minosity seen in both t he MW and M31 dSph 
populations (|Kirbv et al.ll200Sl , l201l"h . 
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Figure 8. Left: r^it vs. cr for the M31 (blue circles), MW (red triangles) and isolated (cyan squares) "classical" dSphs. M31 dSphs 
tend towar ds colder values than the best-fit MW relations (NFW, Cored and power law shown as long-dash, short-dash and solid lines) 
taken from IWalker et all |200§) is observed in ~ 50% of this sample. The "classical" MW dSphs have a tendency to scatter above this 
relation (with the notable exception of CVnl). The two isolated dSphs, Cetus and Tucana are significantly hotter than the MW mass 
profiles. We also show the MW and M31 ultra-faint objects (i.e. My > —7.9) as small green points. The majority of these are consistent 
with the relations, barring And XII and Hercules which are both colder. Right: My vs. r^ a ;/ for the "classical" population, colour-coded 
as before. A number of the M31 dSphs are more extended for a given luminosity than their MW counterparts. Interestingly, a number 
of the more extended M31 objects also fall within the colder regime of the previous plot. Lower: Central surface brightness (S^q) vs - 
scatter for all dwarfs about the best fit NFW profile from the left hand panel, where the scatter is measured as the difference between 
the observed velocity dispersion from the observed velocity dispersion, divided by the expected dispersion in the NFW model, <J v /a n t w . 
If these dSphs all inhabited a "Universal" halo, you would expect to see the points scattering around zero, but instead we see a trend of 
increasingly negative scatter with dec reasing central surface b rightness. We over-plot the expected path for an object undergoing tidal 
disruption to follow in this plane, fron jPenarrubia et al.l (j2008T l . The black points represent a dwarf which has yet to lose any mass from 
tides, and the point at which it has lost nine tenths of its mass. It follows our observed trend well, and we argue that tides have likely 
played some part in the scatter about a "Universal" mass profile. 
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4 COMPARING THE "CLASSICAL" DSPHS 
OF M31 AND THE MW 

In the literature, the term "classical" has been applied to 
the MW dSph galaxies that were discovered pre-2005 i.e. 
Carina (Car), Draco (Dra), Fornax (For), Leo I (Leol), Leo 
II (LeoII), Sculptor (Scl), and Ursa Minor (UMi). All of 
these dSphs have absolute magnitudes of My < —8. By 
limiting ourselves to this brighter regime, we can ensure 
that we have robust estimates for the parameters we are 
interested in, particularly the velocity dispersions. While 
there are a few dSphs in M31 with My > —8 that have 
been spectroscopically surveyed (i.e . And XI, XII and XIII, 
IChapman et all 120071 ; ICollins et aTl 120101 ) . their RGBs are 
sparsely populated, leaving us with only a handful of stars 
(4-8) with which to determine kinematic properties. This 
results in large associated errors. For the "classical" dSphs, 
one is able to obtain velocities for much larger numbers of 
members, giving us velocity dispersions that are much bet- 
ter constrained. Therefore for the purpose of this study, we 
require "classical" dwarfs to be any dSph with My < —7.9. 
With this luminosity constraint, we ensure a quality control 
for our M31 dSphs, and avoid the regime of the ultra-faints 
in the MW where a number of factors (as discussed in the 
introduction) make the determination of the masses and ve- 
locity dispersions in these systems very difficult. Multiple 
epoch data is required for these obj ects in order to con - 
strain their true velocity dispersions l|Koposov et aljfeoi ll ). 
We choose a lower cut of My < — 7.9 rather than My < —8 
as it allows us to include the MW object Canes Venatici I 
(CVnl) in our sample. This object has more in common with 
the classical MW dSphs than the ultra-faint population, and 
it is quite similar to a number of the classical M31 dSphs in 
terms of its luminous component, with a larger rhaij for its 
luminosity than typical MW objects. Therefore it is inter- 
esting to examine how it compares with both populations. 

With the sample d efined, we now investigate the claims 
of IWalker et ail (|2009h that the MW dSphs are consis- 
tent with being born in a "Universal" dark matter halo, 
i.e. a halo where the mass contained within a given ra- 
dius is identical for all dSphs, irrespective of their lumi- 
nosity. We plot rhaif vs. a v , for our sample, adding our 
newly derived velocity dispersions for And V and And VI 
in Fig. [S] where the MW dSphs are plotted as red trian- 
gles, and the M31 objects as blue circles. Here we over- 
lay the best-fit cor ed and Nava r ro-Fre nk- White (NFW) 
mass profiles from IWalker et all (2009), which were de- 
rived using the a v and r haif data from MW satellites. 
Our MW data comes from IWalker et all ([2009), and refer- 
ences within, with an updated value for th e velocity disper- 
sion of Bootes from (Koposov et al.|[201l|). Our M 3 1 dat a 
are taken from iLetarte et al.1 l\200$ ), ICollins et all (|201dh . 
iKalirai et all (|20ld ) and references within. Our two new 
measurements for A nd V and And VI agree well with the 
IWalker et al] l|2009l ') relations, as do And I, And VII and And 
XV, demonstrating that while M31 does appear to have a 
population of colder outliers, among its brighter population 
there is a significant fraction of dSphs that agree well with 
their MW counterparts. This raises the question, what is the 
difference between these systems and their colder brethren? 

Whilst the M31 dSphs scatter towards colder velocity 
dispersion for a given size, the MW dSphs scatter in the op- 



posite direction, with 4 of the 9 ( Dra, Scl, UMi and For) 
satellites lying over la from the NFW and Core profiles in 
the hotter sense, while only one member of this population 
(CVnl) falls below it. Even as we go to lower luminosities, 
shown as small green points, we find only one other dSph 
that is colder than th ese relations, the Hercules dSph, which 
several authors (e.g. ISand et al.ll2009l ; I Jin fc Martini l2010h 
claimed may not be a dSph at all, but a stellar stream of 
debris formed from the disruption of a dwarf satellite. Sim- 
ilarly, while 5 of our 10 M31 objects are significantly colder 
than these relations, only one of these objects (Andl) ap- 
pears to be significantly hotter than predicted. We show the 
two isolated Local Group dSphs, Cetus (Cet) and Tucana 
(Tuc) in Fig. [8] and these are also seen to be outliers to the 
MW relations, with hotter dispersions than expected from 
the best-fit profiles. 

To see how significant these deviations are, we measure 
the statistical offset of the "classical" dSphs of both the MW 
and M31 from these profiles. To do this, we calculate the dif- 
ference between the velocity dispersion of each object, and 
the expected velocity dispersion given by the best-fit NFW 
and Core profiles for a dSph of the same half-light radius, 
and we divide this by the la errors associated with the ob- 
ject in question. This will allow us to measure the scatter 
for both populations, as well as the average offset, (\a\). 
Taking the MW sample first, we measure an average devi- 
ation of (\a M w,NFw\) = 1-64 and (\(TMW,Core\) = 1-44. If 
these profiles were truly a good fit for this population, one 
would expect to see a scatter of order (\a\) ~ 0.2 — 0.3. 
This implies that the MW dSphs in our sample are not well 
fit by these mass profiles. In fact, if we look at the dis- 
tribution of the scatter, we find that 5 out of 9, or 55% 
of the dSphs in this sample are distinct from these pro- 
files beyond their la errors. Further, we note that 2 of the 
MW dSphs, For and CVnl are outliers at > 3<r for both 
profiles. If these profiles were truly Universal, one would 
not expect to see two such significant outliers from a pop- 
ulation of 9 objects (i.e. 22%). We also perform this exer- 
cise in a sign-dependent fashion to see if there is a pre- 
ferred direction for the scatter in the MW dSphs about 
these profiles. Here we find (aMW,NFw) = +0.72 ± 0.2 
and (aMW,Core) = +0.5 ± 0.1, where the errors represent 
the standard deviation of these measurements. This sug- 
gests that the MW sample scatter more in the positive or 
"hotter" direction about these best-fit profiles. Performing 
the same analysis with the M31 "classical" dSphs, we find 
(\aMW,NFw\) = 1-52 and (\aMW,Core\) = 1-4. This shows 
that the of scatter about these profiles for the MW and M31 
dSphs is very similar, and that the profiles of I Walker et al.l 
(2009) are not a demonstrably worse fit to these objects. 
The typical errors on the M31 velocity dispersions are larger 
than their MW counterparts, and as these are reduced in size 
with future observations, we will see if this fit gets better or 
worse. Again, we note that there are 2 significant outliers 
in the M31 system, And II and And X, both of whom sit 
more than 3a below the best fit relations, and that 6 of the 
10 M31 dSphs (60%) are outliers to the relation at > la. 
Again, assessing the scatter of the population in a sign- 
dependent manner, we find {a A/31, nfw) = —0.80 ± 0.2 and 
(o\M3i,Core) = —0.59 ± 0.4 showing that these objects pref- 
erentially scatter in the more negative or "colder" direction. 
Overall, we find that the distribution about these profiles 
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(for both M31 and MW) demonstrates that there is more 
scatter overall than you would expect from solely their indi- 
vidual errors, as only 8 out of 19 (or 42%) objects sit within 
la of their expected values. Finally, we note that the two 
isolated dSphs, Cetus and Tucana, are significant (> 4cr for 
Cetus and > 2a for Tucana) outliers to the relations. 

This spread in the overall population (MW, M31 and 
isolated dSphs combined) suggests that some physical fac- 
tor is driving the scatter about these relations, and we ar- 
gue that tidal forces exerted by the host could play a part in 
creating the observed differences in both the central velocity 
dispersions and surfac e brightnesses between M W and M31 
satellites. The work of lPefiarrubia et al.l (|2010h explored the 
masses and kinematics of the dSph populations for two host 
galaxies with the same total mass, but with stellar discs 
whose masses differed by a factor of two. They found that 
the dSphs orbiting the heavier disc had lower masses and in- 
trinsic velocity dispersions for a given half light radius than 
those in the lower disc-mass system, and that these differ- 
ences were caused by tidal forces exerted from this more 
massive central disc . This mechanism was also discussed in 
IWalker et alj ((2010) as a possible explan ation of the dis- 
crepan cy between the M31 dSphs and the [McGaug h et al.l 
(2007) mean rotation curve. We know from various stud- 
ies (e.g. lHammer et al.l [ 200 7i) that the stellar disc of M31 
is roughly twice as massive as the MW disc, and we do in- 
deed see colder dispersions for a given half-light radius in 
M31 cf. the MW. It is thus plausible that this larger disc-to- 
halo mass ratio has influenced the kinematic properties for 
a number of M31 dSphs. Those that seem more typical in 
terms of their velocity dispersions and central masses could 
then either be on orbits about their host that cause them to 
feel the tidal forces of the disc less keenly than their coun- 
terparts, such as And V and And VI, or could have been 
accreted to M31 at a later time. The tendency for the MW 
objects to spread in the opposite direction could also be ex- 
plained using similar arguments. Tides may have played a 
more significant role in the evolution of the colder CVnl, the 
only classical MW that falls significantly below the Walker 
2009 relations. While no extreme extra-tidal population has 
been reported for this object, it is q uite extended along its 
major axis (an extension of ~ 2kpc. lMartin et alj|2008f ). 

Another discrepancy between brighter MW and M31 
dSphs is that the latter have a tendency to be more ex- 
tended for a given luminosity than their MW counterparts 
(IMcConnachie fe Irwinl [20063). In the right hand panel of 
Fig. [8] we plot the absolute magnitudes of both the "clas- 
sical" M31 and MW dSph as a function of their half-light 
radius. Interestingly, a number of the more extended M31 
objects in this panel also fall within the colder regime of the 
previous plot, meaning the colder objects tend to have larger 
scale radii than those that are observed to be more typical 
or hotter. In the lower panel of Fig. [8] we demonstrate this 
trend by plotting the logarithm of central surface bright- 
ness, log[E v ,o] (where we measure T, v ,o = L ha i f /nrl alf ) 
against the logarithm of the scatter about the best fit NFW 
profile in the a v — r^aif plane, where the scatter is mea- 
sured as the ratio of the observed velocity dispersion to 
the expected dispersion in the NFW model, log[<T v /o"nfw]. 
If these dSphs all inhabited a "Universal" halo, you would 
expect to see the points scattering around zero, but in- 
stead we see a trend of increasingly negative scatter with 



decreasing central surface brightness. We argue that this 
is likely driven in part by tidal forces. As the dark mat- 
ter halos are gradually stripped of mass, they become less 
dense, causing a drop in their velocity dispersions. At the 
same time, the stellar surface density decreases as stars are 
gradually removed. To illustrate the expected path taken by 
dSphs undergoing tidal disruption in this parameter space, 
we over plot the tidal trac k s gene rated from eqn. 7 and Ta- 
ble 2 of iPenarrubia et akl (|2008T ) and note that it follows 
the same direction as our data, indicating that tides have 
likely played some part in the evolution of these objects 
away from their initial mass profiles. The black points rep- 
resent a dwarf which has yet to lose any mass from tides, 
and the point at which it has lost nine tenths of its mass. 
This figure clearly demonstrates that for some of the highest 
surface brightness MW objects (e . g. Le o I, Scl and For), the 
NFW profile from IWalker et all (2009) agrees poorly with 
the observed dispersions for these objects, underestimating 
them by ~ 0.1 dex. These objects would have their central 
dispersions and masses modeled better with halo profiles 
that possessed a hi gher circular velocity than that used by 
IWalker et ail (|2009D (V c ~ 18kms _1 c/.V c ~ lSkms" 1 ), al- 
though such a profile would be demonstrably worse for the 
colder MW and M31 objects. This finding strongly argues 
against the notion that all dSphs currently reside in a "Uni- 
versal" halo. We also plot the positions of the ultra- faint 
dSphs in this plane, however owing to the large uncertain- 
ties in measured values of luminosity, half-light radius and 
velocity dispersion, they demonstrate a large scatter about 
the expected value of zero. We note that there are a num- 
ber of ultra-faints that have significantly colder dispersions 
than expected for their surface brightness, including Her- 
cules, Bootes and And XII, suggesting that tides likely play 
an important role in shaping these objects. 

Exploring this tidal origin as an explanation for these 
differences we turn to the kinematic properties of the two 
isolated Local Group dSphs, Cetus and Tucana. Owing to 
their large distances from either the MW or M31 (755 kpc 
and 890 kpc), they are not thought to have felt a strong 
tidal force from either galaxy over the course of their evo- 
lution. For the Cetus dSph, analysis of both INT photome- 
try an d Keck DEIMOS sp ectroscopy l|McConnachie fe Irwinl 
l2006al ; H ewis et al. 2007) has shown no evidence for previ- 
ous or o ngoing tidal disruption , and it s large tidal radius 
(6.6 kpc IMcConnachie fe Irwinl [2006a) implies that tidal 
forces from either the MW or M31 are unlikely to have 
played a significant role in its evolution. Similarly, Tucana 
shows no obvious signs of tidal disrupti on, other than an 
absence of HI gas (|Fraternali et al.ll2009T ) which could have 
been expelled by stellar feedback, or tidally removed. As 
it has an unusually large rec eding velocity with res pect to 
the MW (v r = +98.9kms- 1 , lFraternan et al.|[2009h . it has 
been argued that Tucana may be on a highly elongated orbit 
about the MW, bringing it in close proximity to our Galaxy 
~ lOGyrs ago, meaning it may have experienced stronger 
tidal forces than Cetus. Given its current position, it is un- 
likely to have completed more than one orbit within the Lo- 
cal Group, and therefore has experienced less tidal-stirring 
over the course of its history compared with the M31 and 
MW dSphs. In our analysis here, we find both these objects 
are positioned in the hotter r^aif — &v regime, with Cetus 
being a more extreme outlier than Tucana. This fits with 
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our expectation if tidal forces from the host are the cause of 
this scatter about a Universal profile. 

These results demonstrate that, while other factors may 
also be involved in shaping the physical and dynamical prop- 
erties of the dark matter halos of dSphs, the effects of tidal 
forces exerted by the host galaxy play an important role in 
setting the underlying mass profiles of its associated satel- 
lite population. In the context of the MW and M31 systems, 
we argue that the observed scatter towards colder velocity 
dispersions in the M31 subhalos is likely due to the larger 
disc-to-halo mass ratio in this system compared with the 
MW. 



5 CONCLUSIONS 

In this study, we set ou t to fu rther investiga t e the results 
of ICollins et all (|2010h and iKalirai et al] |2010l ) who 
noted that the velocity dispersions and central masses of 
a number of the M31 dSphs differed significantly from 
MW dSphs of similar scale-radii, making them outliers 
to the uni versal mass profi l es der ived for the MW dSphs, 



and the McGaug h et al 



(|2007l ) mean rotation curve 
|Walker et alj 120091 . l201fj ). We restricted ourselves to 
studying the "classical" dSphs in both galaxies, where 
classical was defined as any dSph with My < — 7.9 to 
avoid complications introduced by the uncertain nature 
of some ultra-faint dwarfs, and low number statistics 
for the faint M31 dSphs. We also presented updated 
structural properties and kinematic properties for two 
of the classical M31 population, And V and And VI, 
for whom only systemic velocities have previously been 
reported ()Guhathakurta et al.l l200d ). We updated values 
for their half-light and tidal radii to Thai f = 292 ± 22pc and 
r half = 440±16pc, n = 1.2±0.2kpc and r t = 1.6±0.2kpc, 
their PA (0 = 32 ± 2° and 9 = 164 ± 2°) and ellipticities 
(e = 0.17 ± 0.02 and e = 0.39 ± 0.02) for And V and 
And VI respectively. In terms of their kinematics, we mea- 
sured systemic velocities of v r — —393.1 ± 4.2 kms -1 and 
344.8 ± 2.5 kms -1 , and dispersions of a v — II.5I44 km s _1 
and a v = 9. i^'l kms -1 for And V and And VI respectively, 
meaning that these two dSph appear very typical for objects 
of their size when compared with MW dSphs. This result 
shows that not all of the M31 dSphs reside in significantly 
different dark matter halos to those of the MW. When 
assessing the classical MW and M31 in the rhulf — &v plane, 
we find that with r espect to the best-fit mass profiles of 
IWalker et all l|2009h . scatter about these is observed in 
both the positive and negative directions, and is greater 
than would be expected from the measurement errors alone. 
The scatter in the positive or "hotter" direction occurs 
predominantly within the MW population with only one 
M31 object (Andl) found to scatter in this direction, and 
similarly, all but one of the "colder" objects are M31 dSphs 
(the exception being CVnl). We also find a tendency for the 
"hotter" dSphs to be more compact for a given luminosity 
than the "colder" objects. This is seen when comparing the 
central surface brightnesses of these objects w ith deviations 
of the ir velocity dispersions from the best-fit Walker et al.l 
(2009) NFW profile. In this plane, we can clearly see that 
deviations towards colder velocity dispersions increase as 
surface brightness decreases. Analysing this in the frame- 



work of IPefiarrubia et al.l |2010l ) where the colder velocity 
dispersions of the M31 dSphs are suggested to be a result 
of a more massive disc-to-halo mass ratio in M31 compared 
to the MW, we argue that the underlying mass profiles for 
dwarf galaxies are not Universal, and are influenced the 
baryonic component of disc galaxies, causing variation from 
host to host. 
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